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Abstract 


The metastable helium line at 1083 nm can be used to probe the extended upper atmospheres of close-in exoplanets 
and thus provide insight into their atmospheric mass loss, which is likely to be significant in sculpting their 
population. We used an ultra-narrow band filter centered on this line to observe two transits of the low-density gas 
giant HAT-P-18b, using the 200" Hale Telescope at Palomar Observatory, and report the detection of its extended 
upper atmosphere. We constrain the excess absorption to be 0.46% + 0.12% in our 0.635 nm bandpass, exceeding 
the transit depth from the Transiting Exoplanet Survey Satellite (TESS) by 3.9c. If we fit this signal with a 1D 
Parker wind model, we find that it corresponds to an atmospheric mass loss rate between 8.3* a x 10-5 M; Gyr! 
and 2.637046 x 10-3 M; Gyr ' for thermosphere temperatures ranging from 4000 K to 13,000 K, respectively. 
With a J aami of 10.8, this is the faintest system for which such a measurement has been made to date, 
demonstrating the effectiveness of this approach for surveying mass loss on a diverse sample of close-in gas giant 
planets. 


Unified Astronomy Thesaurus concepts: Exoplanet atmospheres (487); Narrow band photometry (1088); Transits 


(1711); Exoplanets (498); Extrasolar gaseous giant planets (509); Extrasolar gaseous planets (2172) 


Supporting material: data behind figure 


1. Introduction 


Close-in exoplanets are exposed to high-energy radiation from 
their host stars, which can lead to atmospheric mass loss. This 
atmospheric escape appears to shape the observed short-period 
exoplanet population (e.g., Lopez & Fortney 2013; Owen & 
Wu 2013; Fulton et al. 2017), but there are relatively few 
published measurements of present-day mass loss rates for close- 
in planets. Prior to 2018, most studies of atmospheric escape used 
the hydrogen Lya line at UV wavelengths, the Ha line at optical 
wavelengths, and metal lines at UV and optical wavelengths (e.g., 
Vidal-Madjar et al. 2003, 2004; Jensen et al. 2012; Yan & 
Henning 2018; Cauley et al. 2019). In recent years, new 
theoretical and observational work on the helium (He) 1083 nm 
line have shown that it can also be used for atmospheric mass loss 
studies. For planets with a sufficient population of metastable 
helium atoms in their (potentially escaping) upper atmospheres, 
the 1083nm line is optically thick at low pressures (high 
altitudes), increasing their measured transit depths in this line by a 
factor of a few (Okloptié & Hirata 2018). K-type stars are 
favorable targets for metastable helium observations, as they emit 
relatively low amounts of mid-UV flux (which depopulates the 
metastable state), while emitting relatively high levels of EUV 
flux (which populates the metastable state via ground-state 
ionization and subsequent recombination; Oklopéié 2019). 

Excess He I absorption was first detected in the atmosphere of 
the sub-Satum WASP-107b using the Hubble Space Telescope 
(HST; Spake et al. 2018). Since then, this line has been used to 
detect extended atmospheres in six other planets (five gas giants 
and one sub-Neptune, with masses ranging from 0.044Mj 
to 1.116My) using both space- and ground-based facilities (Allart 
et al. 2018; Mansfield et al. 2018; Nortmann et al. 2018; 


Salz et al. 2018; Alonso-Floriano et al. 2019; Ninan et al. 2019; 
Palle et al. 2020). Understanding how the metastable helium 
signal scales with stellar EUV flux, planetary gravity, and 
stellocentric distance on a sample of large, well-characterized 
planets is a prerequisite for He I observations of smaller planets, 
where measurements and interpretation are inherently more 
difficult (Kasper et al. 2020). Additionally, detections of 
outflowing gas giant atmospheres in the He I line provide useful 
constraints on crucial radiative and collisional processes in 
theoretical atmospheric mass loss models (e.g., Salz et al. 2016; 
Oklopéié & Hirata 2018). 

Using the He 1083 nm line, we study the extended atmosphere 
of HAT-P-18b (Hartman et al. 2011), which is a Jupiter-sized 
(0.947 + 0.044R;), Saturn-mass (0.196 + 0.008M)), Teg = 841 + 
15 K planet (Esposito et al. 2014) orbiting a K2-type star with 
J=10.8 (Cutri et al. 2003). The low density of the planet 
combined with the spectral type of the host star makes it an 
excellent target for metastable helium studies. Previous Rossiter— 
McLaughlin studies of this system have shown it to be one of the 
only planets around a cool main-sequence star with a retrograde 
orbit (Esposito et al. 2014). Additionally, optical transmission 
spectroscopy has revealed a Rayleigh scattering slope for this 
planet, potentially due to a high-altitude haze (Kirk et al. 2017). A 
secondary eclipse of HAT-P-18b has also been detected by the 
Spitzer Space Telescope, and the resulting brightness temperature 
for the planet suggests efficient day-night circulation and/or a 
nonzero albedo (Wallack et al. 2019). 

In this work, we characterize the atmospheric mass loss of 
HAT-P-18b for the first time. We observed two transits of HAT- 
P-18b with the Hale 200” Telescope at Palomar Observatory 
using an ultra-narrow band filter centered on the helium 1083 nm 
line (Vissapragada et al. 2020) with the Wide-field InfraRed 
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Camera (WIRC; Wilson et al. 2003). Additionally, we use TESS 
data from Sectors 25 and 26 to improve the ephemeris for this 
planet and determine its broadband optical transit depth as a 
comparison for our helium measurement. In Section 2 we describe 
the WIRC and TESS observations, and in Section 3 we jointly 
model the light curves from both instruments. We describe our 
results in Section 4, and offer some concluding thoughts in 
Section 5. 


2. Observations 
2.1. WIRC Observations 


We observed transits of HAT-P-18b on UT 2020 June 6 
(hereafter night 1) and 2020 July 8 (night 2). We used an ultra- 
narrow band filter centered on the helium 1083.3 nm line with 
an FWHM of 0.635 nm (Vissapragada et al. 2020). For these 
observations we typically use a custom beam-shaping diffuser, 
which produces a stable 3" diameter top hat point-spread 
function (Stefansson et al. 2017), to mitigate time-correlated 
systematics (Vissapragada et al. 2020). However, on night 1 the 
weather conditions were poor and we elected to defocus the 
telescope to 1/2 rather than using our diffuser, as the precision 
of the photometry would have been compromised by the 
increased sky background introduced into the diffused point- 
spread function of the star. On night 2, we utilized the diffuser 
as normal. Before each observation we observed a helium arc 
lamp installed at the Hale 200" to determine our region of 
maximal sensitivity to the 1083 nm line, and we placed our 
target on this region during data collection. Due to the change 
in the filter's center wavelength across the detector, telluric OH 
emission lines form bright radial arcs on the detector; to 
calibrate this we constructed a background template for the OH 
lines using a four-point dither on each night. For our night 1 
observations, we took 90 s exposures from UT 05:03:30 to UT 
10:40:11, beginning at airmass 1.237 and ending at airmass 
1.207. For our night 2 observations, we took 90 s exposures 
from UT 05:06:57 to UT 11:01:56, beginning at airmass 1.011 
and ending at airmass 2.301. Both nights reached a minimum 
airmass of 1.000. 

To calibrate the images, we dark-subtracted and flat-fielded 
all the science data, correcting for bad pixels and any residual 
detector striping. This process is described in Vissapragada 
et al. (2020). To correct for the bright arcs on the detector 
caused by telluric OH, we median-scaled the sigma-clipped 
science data to the dithered background frame in 10 pixel radial 
steps from the filter zero-point (where rays have passed through 
the filter at normal incidence) at the top of the detector. This 
process eliminates most of the telluric background, leaving a 
small residual background that we correct for locally in our 
aperture photometry process. 

We use these background scaling factors to correct for time- 
varying telluric water absorption as well. There are two water lines 
that overlap with the bandpass of our filter at 1083.57 and 
1083.66 nm. At the effective resolving power of our filter, these 
two lines appear as a singular absorption line. Previously in 
Vissapragada et al. (2020), we could assume that the strength of 
the absorption in this water line was effectively constant over the 
night, as the observing conditions were good and the resulting 
light curves did not exhibit time-correlated variations corresp- 
onding to a rapidly varying water vapor column. However, since 
the first of our nights had relatively poor weather conditions with 
sporadic cloud coverage and noticeable seeing and transparency 


Paragas et al. 


variations, we sought to track variations in this water absorption 
feature using the OH sky emission lines. The Lorentzian wing of 
the Q,(3/2) OH emission line at 1083.4 nm overlaps with the 
telluric water feature (Allart et al. 2018; Salz et al. 2018). This 
means that OH emission originating higher up in the atmosphere 
(280 km; Bernath & Colin 2009) can be absorbed by H2O while 
passing through the lower atmosphere. We can therefore track the 
telluric water variation over the night by dividing the time-varying 
flux in the water-contaminated OH emission line (as measured by 
our scaling factors) by that of the uncontaminated R,(3/2) and 
R»(1/2) OH lines at 1075 and 1078 nm. If the water variation is 
significant enough to impact our photometry, we can utilize this 
absorption proxy as a decorrelation parameter in our transit fits. 

We performed aperture photometry on the target star and six 
comparison stars (the same ones on each night) using the 
package photutils (Bradley et al. 2016). We tested 
different aperture sizes in one pixel steps from 3 to 13 pixels 
in radius. We removed 3c outliers from the data using a 
moving median filter. This process is detailed in Vissapragada 
et al. (2020). Our optimal aperture sizes for our night 1 and 
night 2 observations were 7 pixels (1"75) and 11 pixels (2775) 
in radius, respectively, with the difference arising from our use 
of a slight defocus with the bandpass of our filter on the first 
night and a diffuser on the second night. The raw light curves 
of the target and the comparison stars are shown in Figure 1 for 
both nights. 


2.2. TESS Observations 


We used the 2 minute cadence TESS observations of HAT- 
P-18b obtained during Sectors 25 and 26. TESS observed the 
target for 51.5 days starting on 2020 May 14 and 2020 June 9 
for Sectors 25 and 26, respectively, covering eight transits in 
total. We downloaded the Pre-search Data Conditioning Simple 
Aperture Photometry (PDCSAP) light curve from the Mikulski 
Archive for Space Telescopes (MAST) using the light- 
kurve package (Lightkurve Collaboration et al. 2018). With 
the transits masked, we removed low-frequency variability 
from the data using the Savitzky-Golay filter from scipy 
(Virtanen et al. 2020) and rejected 5c outliers using a moving 
median filter. However, we noticed that even after the filter was 
applied, there were still strong uncorrected systematics that 
biased the transit depths of the first two transits, so we omitted 
them from our combined fit. Although these transits may be 
recoverable with different detrending methods, our constraint 
on the TESS transit depth from the remaining six transits were 
sufficiently precise for comparison to the WIRC light curves 
(i.e., the uncertainty in the comparison is dominated by the 
uncertainty on the WIRC transit depth). 


3. Light-curve Modeling 


We simultaneously fit both nights of WIRC data along with the 
corrected TESS photometry using exoplanet (Foreman- 
Mackey et al. 2020). For each WIRC light curve, we fit an 
instrumental noise model consisting of a linear baseline along 
with a linear combination of comparison star light-curves, with the 
weights of the comparison stars left as free parameters in the fit. 
This is an update from our previous modeling methodology 
(Vissapragada et al. 2020), where we used an ordinary least- 
squares method to quickly determine comparison star weights at 
each likelihood evaluation, and it is enabled by the rapid No 
U-Turn Sampler (NUTS) sampler that exoplanet makes 
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Figure 1. (a) Results for WIRC night 1 and night 2 shown in (a) and (b), respectively. The top row shows the median-normalized raw light curves of the target (in 
black) and its comparison stars (in gray) as well as the water absorption proxy described in Section 2.1 (dotted blue line) throughout each respective night. The middle 
row shows the helium light curve with unbinned data in gray and binned data to a 15 minute cadence in black, with the best-fit joint helium model in red and the TESS 
model in blue, and the residuals of each fit. The bottom row shows the Allan deviation plot for each data set. 


available for high-dimensional light-curve fitting. We also tried 
including two additional parameters in our instrumental noise 
model for each night: the water absorption proxy (as described in 
Section 2.1), and the distance from the median centroid. We find 
that two comparison stars in each WIRC night have posterior 
probability distributions for their weights that overlap with zero, 
and we therefore remove them from the fit, lowering the Bayesian 
information criterion (BIC) value by 35 and 28 for night 1 and 


night 2, respectively. Although these stars are not the same for 
each data set, the two nights had different observing strategies and 
weather conditions. We find that we obtain optimal fits when we 
only include the telluric water proxy in our fits to the night 1 data, 
lowering the BIC value by 17; this is not surprising, as this night 
had relatively poor and variable weather conditions. We opted to 
keep the distance from the median centroid as a decorrelation 
parameter for both nights, as their removal from the fits resulted in 
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Figure 2. Combined transit light curves and residuals for WIRC (left) and TESS (right), with unbinned data in gray and binned data to a 15 minute cadence in black. 
The best-fit models for WIRC (red) and TESS (blue) are overplotted with the 1o confidence interval denoted by the shaded region. 


(The data used to create this figure are available.) 


Table 1 

Priors and Posteriors for Joint Fit to Palomar/WIRC and TESS Data 
Parameter Prior Posterior Units 
(Rj/R,? (He) uU, 25) 2.297012 % 
(R,/R,)? (TESS) U, 25) 1.8327 9.08 % 
P N(5.5080291, 0.0000042) 5.508029 + 0.0000042 days 
To (2038.5, 2039.0) 2038.82530 + 0.00023 BTJDrpp 
b N(0.352, 0.057) 0.338*0027 e 
u; (He) Kipping (2013) 0.587020 
u (He) Kipping (2013) 0.14 + 0.39 
u, (TESS) Kipping (2013) 0.45 + 0.16 
u (TESS) Kipping (2013) 0.201939 
log(Gextra) (night 1) u(—4, —2) —2.078* 0036 
log(extra) (night 2) u(—4, —2) —2.407015 
k (TESS) (0.5, 1.5) 0.8563* 0.003 
absorption proxy (night 1) N(0.0, 0.1) —0.078*0026 


Note. BTJD1ps = BJD—2457000. Note that we omitted the stellar parameters and all of the detrending weights except for the absorption proxy for night 1. 


a ABIC< 10. Our final systematics model contained 15 
parameters: two parameters for each of the linear baselines, four 
comparison stars for each data set, the distance from median 
centroid for each data set, and the absorption proxy for the first 
night. 

We fit a transit model simultaneously with the systematics. We 
have three fit parameters that are common to all data sets: the 
predicted mid-transit time 75, the period P, and the impact 
parameter b. Initially, we allowed each night of WIRC data to 
have its own transit depth in the joint fit. The two transit depths 
(Q.11*02396 for night 1 and 2.35 +0.14% for night 2) were 
within 1c of each other, indicating that the magnitude of helium 
absorption appears consistent between these two epochs. We 
therefore fit a single transit depth for both helium light curves. We 
fit for the limb darkening coefficients [uj, U2](He) and [ui, U2] TEss) 
and transit depths (R, [R Juo) and (R, /R eres) in each 
bandpass, both sampled uniformly (see Table 1). For each WIRC 


light curve, we fit for a jitter parameter describing the excess noise 
in addition to the photon noise log(Cextra). For TESS, we noticed 
that the error bars that came with the PDCSAP fluxes were not an 
accurate representation of the photon noise and in fact over- 
estimated the observed scatter in the data. We therefore include a 
scaling factor k for the TESS error bars. 

We use the NUTS in PyMC3 (Salvatier et al. 2016) to sample 
the posterior distributions for our model parameters. We ran 
four chains, tuning each for 1500 steps (the “burn-in” period) 
and then taking 1000 draws in each chain, achieving good 
convergence with a Gelman-Rubin (Gelman & Rubin 1992) 
statistic of « 1.006 for all parameters. The priors and posteriors 
for the physical parameters in our model are given in Table 1 
for the joint fit and the detrended light curve, residuals, and 
Allan deviation plot for each night of WIRC data are displayed 
in Figure 1. The final combined helium and TESS light curves 
are displayed in Figure 2, and the posterior distributions for the 
model parameters are visualized in Figure 3. 
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Figure 3. Corner plot displaying the posterior probability distributions for the joint model for HAT-P-18b. Note transit depth (R,,/ R,,)° values are in %, period P is in 
days, and predicted mid-transit time To is in BTJDrpg. We omit all of the decorrelation parameters except for the absorption proxy. 


4. Results and Discussion 


We measure a transit depth of 2.291012 % in the helium line. 


The corner plot for our fit parameters is shown in Figure 3. Our 
measurement can be compared to the TESS transit depth 
measurement of 1.83270. %, Our measured transit depth in 
the helium bandpass exceeds that in the TESS bandpass by 
0.46% + 0.12% (3.90). The TESS bandpass is between 600 nm 
to 1000 nm, making it reasonable to use as a comparison for 
our measurement. Previous studies on our target have shown 
variation within this range to be limited to variations on the 
order of the scale height (Kirk et al. 2017). The difference 


between the two transit depths exceeds, by an order of 
magnitude, the expected change in transit depth for a one 
(lower-atmospheric) scale-height change in planet radius 
(0.03%) for this target. Thus, the observed excess absorption 
cannot be explained by broad absorption features—for 
instance, by water—in the lower atmosphere. The helium line 
is near an opacity minimum of water anyway, so this 
explanation is disfavored a priori. We conclude that the 
observed absorption indeed arises from metastable helium in 
HAT-P-18b’s extended atmosphere. 

We use the model described in Oklopéié & Hirata (2018) to 
convert our measured excess absorption into a joint constraint 
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Figure 4. Atmospheric mass loss model for HAT-P-18b. Each point is a 
different mass loss model corresponding to specific Ty and M values, and the 
shading indicates the compatibility between the model and our observed excess 
absorption (with the lighter regions indicating the most concordant models). 


on HAT-P-18b's mass loss rate M and upper atmospheric 
temperature Ty. This model calculates the velocity and density 
profiles of a 9096/1090 H/He 1D Parker wind as a function of 
M and To, and then calculates the level populations for helium 
given a UV stellar spectrum. We use the MUSCLES UV 
spectrum of e Eridani (France et al. 2016; Loyd et al. 2016; 
Youngblood et al. 2016), which is another K2 type star, as a 
stand-in for the unknown UV spectrum of HAT-P-18. 


Accounting for the stellar radius and semimajor axis of HAT- 


P-18b, the EUV irradiance of. the planet was 8Wm ^ 
integrated between 5.5 and 911 A. The results are shown in 
Figure 4. HAT-P-18b's mass loss rate is likely between 
8.313 x 1075 and 2.637049 x 10 M; Gyr ! for thermo- 
sphere temperatures between 4000 and 13000 K, respectively. 
Using the EUV irradiance above along with an efficiency 
parameter €= 0.1, we can also calculate an energy-limited 
mass loss rate for HAT-P-18b (e.g., Murray-Clay et al. 2009): 


. o ERRPF. 
M = — 28 a4 x 109 g s-l. (1) 
GM, 


This estimate agrees well with the inferred mass loss rates in 
Figure 4 (with uncertainties of a factor of a few accommodated 
by similar uncertainties on the efficiency parameter), suggest- 
ing that our observationally derived constraints are energeti- 
cally feasible. Because e Eridani is a relatively young, active 
star, with log(Ryx) = —4.51 compared to log(Ryx) = —4.80 
for HAT-P-18, we repeated the modeling using the MUSCLES 
spectrum of HD 40307, a fairly inactive (log(Ryx) = —4.99) 
K2.5V star. By trying proxy stars with activity levels on either 
side of HAT-P-18b’s activity, we can get a sense for the 
uncertainty on the result based on our choice of EUV proxy. 
With HD 40307 as a proxy, we found a best-fit 
log(M) = 9.60*01) at 4000 K, and 11.20'017 at 13000K, 
nearly identical to our findings for e Eri. This is because the 
mid-UV to EUV flux ratios between the two stars are quite 
similar (Oklopéié 2019). Although e Eri is a much stronger 
X-ray emitter, the cross section to X-ray photoionization of 
helium is very small compared to the cross section in the EUV 
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near the 504A threshold, so the contribution to the flux- 
averaged photoionization cross section is negligible (Oklopčić 
& Hirata 2018). 

We note that there is a strong degeneracy between the mass 
loss rate and thermosphere temperature due to the complex 
dependence of the outflow velocity and density on the 
temperature and the mass loss rate. This degeneracy could be 
partially resolved with a precise line shape measurement, but 
we do not resolve the line shape in these observations. Due to 
the faintness of HAT-P-18, spectrographs on all but the largest 
telescopes may have difficulty resolving the line shape 
precisely enough to break the degeneracy. Additionally, our 
helium light curve is symmetric across our best-fit mid-transit 
time. However, we cannot exclude the possibility of an 
extended egress, as our combined light curve lacks the 
precision required to significantly detect a trailing helium tail 
for such a faint target. 


5. Conclusions 


In this work, we use an ultra-narrow band helium filter 
centered on the 1083 nm line to observe two transits of HAT-P- 
18b. We detect 0.46% + 0.12% excess helium absorption in the 
planet’s upper atmosphere. This detection corresponds to an 
atmospheric mass loss rate between 8.3'T§ x 1075 and 
2.63046 x 10-3 Mj; Gyr, which means HAT-P-18b is losing 
less than 296 of its mass per Gyr. This is typical for close-in gas 
giants, with other helium outflow detections having mass loss 
rates less than 596 per Gyr (Allart et al. 2018; Mansfield et al. 
2018; Spake et al. 2018; Alonso-Floriano et al. 2019). 

Of the handful of planets with detected helium outflows, 
WASP-107b is the most comparable to HAT-P-18b with a 
similar radius of 0.94 Ry, mass of 0.12 Mj, separation of 
0.55 au, and equilibrium temperature of 770 K (Anderson et al. 
2017). If we assume that HAT-P-18b has a He line shape 
similar to that of WASP-107b, we can invert our excess helium 
transit depth to obtain an estimate of the underlying predicted 
line depth of 4.5% + 1.3%. This is noticeably smaller than the 
7.26% + 0.24% depth measured by CARMENES and Keck/ 
NIRSPEC for WASP-107b (Allart et al. 2019; Kirk et al. 
2020). The difference may be due to the smaller gravitational 
potential of WASP-107b (the mass of WASP-107b has recently 
been suggested to be even lower by Piaulet et al. 2021), or 
differences in the EUV spectra of the two stars (WASP-107 is a 
K6 star while HAT-P-18 is a K2). Detailed comparative 
modeling of these two planets may make clear the primary 
control on the metastable helium signal. 

This is the faintest system (J = 10.8) with detected helium 
absorption thus far, establishing the effectiveness of our 
technique for observing such targets with a mid-sized 
telescope. For reference, the next faintest system with detected 
excess helium absorption is WASP-107b with a J magnitude of 
9.4. Of the 11 planets identified in Kirk et al. (2020) as 
promising targets for observations of helium outflows, many 
are challenging targets with J> 10.5. Our photometric 
technique allows us to begin surveying planets around such 
faint stars, expanding the sample of planets with measured 
metastable helium absorption. Further population-level studies 
of extended atmospheres in the He 1083 line will greatly 
improve our ability to calibrate the mass loss models used to 
elucidate the long-term evolution of the close-in exoplanet 
population. 
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